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Simple curvaton models can generate a mixture of of correlated primordial adiabatic and isocur- 
vature perturbations. The baryon and cold dark matter isocurvature modes differ only by an obser- 
vationally null mode in which the two perturbations almost exactly compensate, and therefore have 
proportional effects at linear order. We discuss the CMB anisotropy in general mixed models, and 
give a simple approximate analytic result for the large scale CMB anisotropy. Working numerically 
we use the latest WMAP observations and a variety of other data to constrain the curvaton model. 
We find that models with an isocurvature contribution are not favored relative to simple purely 
adiabatic models. However a significant primordial totally correlated baryon isocurvature perturba- 
tion is not ruled out. Certain classes of curvaton model are thereby ruled out, other classes predict 
enough non-Gaussianity to be detectable by the Planck satellite. In the appendices we review the 
relevant equations in the covariant formulation and give series solutions for the radiation dominated 
era. 



I. INTRODUCTION 



Recent detailed measurements of the acoustic peaks in 
CMB anisotropy power spectrum by the WMAP satel- 
hte P, are consistent with the standard model of a 
predominantly adiabatic, approximately scale invariant 
primordial power spectrum in a spatially flat Universe. 
Frequently it is assumed the initial power spectrum is en- 
tirely adiabatic, though there is still no compelling jus- 
tification for this assumption. Although adiabatic per- 
turbations are predicted from single field models of in- 
flation 1^ , if one allows the possibility of multiple fields 
in the early Universe then there is also the possibility of 
isocMTOfltere perturbations (also known as entrom 
turbations) |lllllllllll[l3,|ll|ll|ll|llli| 
In particular, the recently proposed curvaton model uses 
a second scalar field (the 'curvaton') to form the pertur- 
bations mill El Us mill. The motivation for this 
is it makes it easier for otherwise satisfactory particle 
physics models of inflation to produce the correct primor- 
dial spectrum of perturbations ■ Various candidates 
for the curvaton have been proposed pi I2I I2II27L . 
A curvaton mechanism has also been considered in the 
pre big bang scenario 2£, .3Q„ ,31^ ,32s, .22.) where it can be 
used to produce an almost scale invariant spectrum. 

The curvaton scenario also has the feature of being able 
to generate isocurvature perturbations of a similar mag- 
nitude to the adiabatic perturbation without fine tuning, 
and therefore is open to observational test. 

Early studies of non-adiabatic perturbations, either 
consider ed p urely isocurvature cold dark matter pertur- 
bations |3J| or mixtures of adiabatic and uncorrelated 
cold dark matter isocurvature perturbations [ssl Is^ Is^ 



18^ . However, as first realized by Langlois 9], the adi- 
abatic and isocurvature components can be correlated 
and this correlation may have interesting observational 
consequences [39]. In Ref. |40j they identified four regu- 
lar isocurvature modes, which in general can have arbi- 
trary correlations with each other and with the adiabatic 
mode. Such general models have many degeneracies and 
are badly constrained by pre- WMAP data ^l!, 42]. De- 
tailed CMB polarization data is expected to help with 
this [43. In Ref. (following Ref. [4| with pre- WMAP 
data) they considered a CDM isocurvature mode with 
an arbitrary correlation to an adiabatic mode and found 
that though not favored by the data, a significant isocur- 
vature contribution was still permitted. Constraints on a 
specific model that doesn't produce isocurvature modes 
were given in Ref. jl^ . 

Here we start in Section ^ by making some general re- 
marks about mixed isocurvature models, and discuss the 
corresponding CMB power spectra predictions. Then in 
Section HTH we discuss current observational constraints 
on totally correlated (or anti-correlated) adiabatic and 
isocurvature perturbations, as predicted by the curva- 
ton model. Various scenarios within the curvaton model 
predict specific ratios of adiabatic and isocurvature per- 
turbations, and can be tested directly. In general we find 
constraints on when the curvaton decayed. 

We use the CMB temperature and temperature- 
polarization cross-correlation anisotropy power spectra 
from the WMAP^ P, ^ observations, as well as 
seven almost independent temperature band powers from 
ACBAR^ 48^ on smaller scales. In addition we use 
data from the 2dF galaxy redshift survey i HST Key 
Project [53, and nucleosynthesis [sj using a slightly 
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modified version of the CosmoMC^ Markov- Chain Monte 
Carlo program, as described in Ref. (52j | . 

For simphcity we assume a flat universe with a cos- 
mological constant, uninteracting cold dark matter, and 
massless neutrinos evolving according to general relativ- 
ity. 



II. PRIMORDIAL PERTURBATIONS AND THE 
CMB ANISOTROPY 



It is well known that the curvature perturbation, in 
the constant density or comoving frame (gauge*) is con- 
served on super-Hubblc scales for adiabatic perturbations 
[53I I54. ISSllSa . ISTj l . This is not the case in the presence of 
isocurvature modes since these source changes to the cur- 
vature perturbation. However, as shown in Ref. [s^ (and 
reviewed in Appendix^), in the presence of isocurvature 
modes the large scale evolution can still be analysed eas- 
ily using the curvature perturbation in the frame in which 
the density is unperturbed, C. This can be expressed in 
terms of the curvature perturbations in the frames in 
which individual species are unperturbed (^t using 



c = 



EM 



(1) 



where the dash denotes the derivative with respect to 
conformal time. For non-interacting conserved particle 
species the individual Ci are conserved on large scales 
if there is a definite equation of state Pi = Pi{pi). In 
this case the evolution of C follows straightforwardly from 
Eq. 1^ depending on the evolution of the background en- 
ergy densities. An adiabatic perturbation is one in which 
Ci = C for *j in which case ^ is constant in time on 
lar ge s cales. The isocurvature perturbations are defined 
as I53I 



5,,.3(C.-0)^-37^(^-^ 



(2) 



where 8p^ = pil^i (no sum) is the density perturbation 
in any frame and Tl is the conformal Hubble rate. We 
consider a fluid consisting of photons (7), massless neu- 
trinos (z/), cold baryons (6) and cold dark matter (CDM, 
c) , where it is conventional to describe the perturbations 
with j = 7, in which case the second index can be omit- 
ted so Sb = Sb,'y, etc. The isocurvature perturbations 
are conserved on large scales where the photon-baryon 
coupling is unimportant. The d are related to the frac- 
tional density perturbations in the unperturbed curva- 
ture frame by = 3{1 + pi/ pi)Ci, and for matter with 
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* In the context of this article the term frame and gauge are effec- 
tively interchangeable. See Appendix^for further discussion. 



constant equation of state are also conserved on large 
scales. 

In general isocurvature perturbations give rise to per- 
turbations in the density, and the universe is no longer 
exactly FRW. One exception to this is when the two mat- 



ter perturbations exactly compensate, so 6p^ 



-Spt 



which case the total matter density is unperturbed, and 
hence the universe evolves as though there were no per- 
turbations. In such a universe the CMB anisotropy would 
be dominated by tiny small scale linear effects due to non- 
zero pressure of the baryons or dark mater, and second 
order effects due to the perturbation in the electron num- 
ber density associated with the baryons. At linear order 
Sp^ = —Spfj, C = Sp^ — Sp^ = is a time independent 
solution to the pressureless perturbation equations, and 
adding this solution to any other solution will make no 
difference to the linear CMB anisotropy or matter power 
spectrum. It follows that an initial isocurvature pertur- 
bation with Sb — l,Sc — is observationally essentially 
indistinguishable from one with Sb = 0,Sc — Pb/Pc- 

We now derive an approximate analytic form for the 
large scale CMB temperature anisotropy in the presence 
of primordial isocurvature and adiabatic perturbations. 
Neglecting a local monopole and dipole contribution, tak- 
ing recombination to be instantaneous at at conformal 
time r*, and assuming the reionization optical depth 
is negligible, the monopole and Integrated Sachs Wolfe 
(ISW) contributions to the temperature anisotropics due 
to scalar perturbations are given by 



ST 



20 



(f)' dX, 



(3) 



where is the Weyl potential (see Appendix IXjl and tq 
is the conformal time today and the integral is along the 
comoving photon line of site. Additional terms which 
arise due to the velocities, quadrupoles and polarization 
at last scattering are generally sub-dominant on large 
scales. Since the pressures are assumed to be zero the 
perturbations are purely adiabatic in the matter era, and 
hence C is constant on large scales. During matter domi- 
nation the potential evolves as = Ci -I- C2/t^, and can 
be related to C using Eq. I)A6|I when the anisotropic stress 
is negligible 



4c + % 



(4) 



Since the radiation to matter density ratio only falls off 
as when the matter dominates, when the approx- 

imation of matter domination is accurate it should also 
be valid to neglect the decaying mode oc l/r^ and as- 
sume (j) « ~fC- Using A^ = 4C^, and neglecting the 
ISW contribution Eq. ||2J) then becomes 



T 5^ 
where from Eq. |^ in matter domination 

C ~ RbCb + RcCc- 



(5) 



(6) 
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Here we define the matter fractions i?b = pb/pm, Rc = 
Pc/Pm, where Pm = Pc + Pb- 

During early radiation domination Q'"^'^ w R-yCj + R^Cu 
from Eq. (Q, where we define the radiation density frac- 
tions Ri, = PuKp-y + Pv),R'y = P-yKP-y + Pv)- Using 

Eqs. 121) and ^ and the constancy of the large scale 
we can therefore relate the large scale temperature 
anisotropy to the primordial adiabatic and isocurvature 
perturbations 

^-^^-Ir'' ~l{RcSc + RbSb) + ^RuS,. (7) 
J 5 5 15 

We can take C'^'^ as a measure of the primordial adia- 
batic perturbation. So this formula shows the effect of a 
mixture of adiabatic and isocurvature perturbations on 
the observed large scale CMB temperature anisotropy. 
This result agrees with that in [s^ , despite errors in their 
derivation which arise from an invalid ansatz for the time 
evolution of the velocities and anisotropic stress (demon- 
strated by counterexample in Appendix Ib)) . However in 
matter domination the velocities are negligible so the er- 
ror is harmless, and for the adiabatic and neutrino modes 
the assumption is correct to the required order during ra- 
diation domination. However unfortunately their general 
result for the evolution of the potential is incorrect and 
cannot be used to improve on the above much simpler 
argument. 

This analytic argument shows the main qualitative 
features, though in reality recombination is far from 
being completely matter dominated, and the ISW and 
other contributions will not be negligible. It is however 
straightforward to compute the CMB and matter power 
spectra numerically [Sfl starting from a series solu- 
tion in the early radiation dominated era (Appendix 0. 

III. CONSTRAINING THE CURVATON MODEL 

The curvaton scenario provides a mechanism for allow- 
ing the inflation potential to have more natural proper- 
ties, at the expense of introducing an additional unidenti- 
fied scalar field which generates the perturbations. In the 
curvaton model the inflaton field drives the initial expan- 
sion and generates an era of radiation domination after it 
decays. The expansion rate then slows and the curvaton 
field can reach the minimum of its potential and start to 
oscillate. During oscillation the curvaton field acts effec- 
tively like a matter component, and its perturbation acts 
like a matter isocurvature mode. As the radiation red- 
shifts further the equation of state then changes to matter 
domination as the curvaton density comes to dominate. 
As the background equation of state changes, a curvature 
perturbation is generated from the isocurvature mode. 
The curvaton then decays into (predominantly) radia- 
tion well before nucleosynthesis, and we enter the usual 
primordial radiation dominated epoch. 

Primordial correlated isocurvature modes can be gen- 
erated if the baryons or CDM are generated by, or be- 



fore, the curvaton decays, as discussed in detail below. 
If one or both were created before the curvaton decays, 
the current model assumes that the curvaton had a neg- 
ligible density when they decayed We assume that 
the curvaton is the only cosmologically relevant scalar 
field after infiaton decay, and that the perturbations in 
the infiaton field are negligible. Generically such mod- 
els predict a very small tensor mode contribution, which 
we assume can be neglected. We assume there is no lep- 
ton number at neutrino decoupling so that there are no 
neutrino isocurvature modes, though see for other 
possibilities. 

As discussed in Section^ the baryon and CDM isocur- 
vature modes predict proportional results, so we can ac- 
count for Sc by using just an effective baryon isocurvature 
perturbation 

Sf^Sb + ^S,. (8) 

The baryon and CDM isocurvature perturbations are 
completely correlated (or anti-correlated) with each other 
and the adiabatic perturbation, so Sf^ — B^^^'^ where 
B measures the isocurvature mode contribution and is 
taken to be scale independent^. From Eq. lO, the large 
scale CMB anisotropy variance is then given approxi- 
mately by 

where Vq — 'Pp'^ is the initial power spectrum. We as- 
sume V( is well parameterized by Vq — As{k/ko)"-'~^ 
where As gives the normalization, is the scalar spec- 
tral index and kg is a choice of normalization point. Note 
that our number of degrees of freedom is actually less 
than generic infiation, because although we have intro- 
duced B we now no longer have the amplitude and slope 
of the tensor component to consider. The slope of the 
isocurvature perturbation is predicted to be the same as 
the adiabatic perturbation and the tensors are predicted 
to be negligible in the curvaton scenario [T^ . 

The isocurvature modes have little effect on small 
scales, but as can be seen from Eq. they can either 
raise or lower the Sachs Wolfe plateau relative to the 
acoustic peaks depending on the sign of B. This is in 
contrast to tensor perturbations which can only raise the 
Sachs Wolfe plateau relative to the acoustic peaks. 

Computing the full predictions numerically and assum- 
ing a fiat prior on B, Fig. shows the posterior dis- 
tribution for the various cosmological parameters when 
the possibility of a totally correlated mixture of matter 



^ Our sign convention for B differs from that in Ref. l44l . In our 
convention B > corresponds to a positive correlation and the 
modes contribute with the same sign to the large scale CMB 
anisotropy. 
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FIG. 1: Posterior marginalized probability distributions 
(solid lines) of the cosmological parameters including corre- 
lated matter isocurvature modes, using the data described in 
the text. B is the ratio of the (effective) baryon isocurvature 
to adiabatic perturbation amplitude in the primordial era, 
Qbh^ and Qch^ are the physical matter densities in baryons 
and CDM, Ho km s~^Mpc~^ is the Hubble parameter today, 
Zie is the effective reionization redshift, and Us is the spectral 
index. We assume a flat universe with cosmological constant. 
Dotted lines show the mean likelihoods of the samples, and 
agree well with the marginalized curves, indicating the full 
distribution is fairly Gaussian and unskewed 62]. 
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FIG. 2: Posterior distribution of B = <SS; /(^ in the primordial 
era, and the spectral index Us- The plot is generated from a 
smoothed number density of Monte Carlo samples generated 
using the data described in the text. The contours enclose 
68% and 95% of the probability, and the shading is by the 
mean likelihood of the samples. 



isocurvature and adiabatic perturbations is allowed. The 
posterior distribution of B and is shown in Fig. |21 
marginalized over the other parameters. On small scales 
the isocurvature modes have only a small effect, so the 
main observational constraint comes from the relative 
amplitudes of the large and small scale power. This is 
partially degenerate with the spectral index as clearly 
demonstrated in the figure. The relative large scale am- 
plitude is also affected by the reionization optical depth, 
and although this is constrained by WMAP's polariza- 
tion measurements the experimental noise and cosmic 
variance still leave a significant residual uncertainty. 

We find the ratio of the mean likelihood allowing for 
isocurvature modes to that for purely adiabatic mod- 
els is about 0.7 (for discussion of mean likelihoods see 
Ref. |52]). Thus the isocurvature modes do not improve 
the already good fit to the data of the standard purely 
adiabatic case. By the same token, the current data 
is still consistent with a significant isocurvature contri- 
bution, with the 95% marginalized confidence interval 
-0.53 < B < 0.43. If new data favored B > this 
would be largely degenerate with a tensor contribution 
predicted by standard single field inflationary scenarios, 
and would be hard to distinguish without good CMB 
polarization data. Evidence for i? < would be a smok- 
ing gun for an isocurvature mode, though the large scale 
polarization data has large enough cosmic variance that 
to distinguish it from an adiabatic model with an unex- 
pected initial power spectrum shape would be difficult. 

The 95% confidence marginalized constraint on the 
spectral index 0.90 < Ug < 1.02 translates into a con- 
straint on the potential V during inflation (in general a 
function of the inflation field 'ijj and the curvaton field 

a) (a 

-0.1<2(r;,, -e)<0.02 (10) 

where 

Mp is the reduced Planck mass, and the quantities are 
evaluated at horizon crossing during inflation. In stan- 
dard inflationary models the potential has to satisfy 
1/1/4 ^ 0.03e^/^Mp to obtain the correct fluctuation am- 
plitude, which is difficult without using unnatural values 
of the model parameters |23|. In the curvaton scenario 
we assume the inflaton perturbations are negligible, and 
hence the potential merely has to be much smaller than 
this number. These conditions are therefore much easier 
to satisfy with natural values for the model parameters 
in the curvaton case . In both cases the inflaton com- 
ponent of the potential also has to provide more than 
about 60 e-folds of inflation. 

If the CDM is created before the curvaton decays, and 
while the curvaton still has negligible energy density, its 
density is essentially unperturbed. After the curvature 
perturbation is generated there is therefore a relative 
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isocurvature perturbation, given by |19| 
Sc « -3C. 



(12) 



If the curvaton decays before its energy density com- 
pletely dominates, a CDM isocurvature perturbation is 
produced 



1 



(13) 



where r measures the transfer function from Ccurvaton 
before curvaton decay to C after decay, C = ''Ccurvaton- 
Ref. ^3 find the approximate result r « Pcurvaton/ptotai 
where p is the energy density at curvaton decay, to an ac- 
curacy of about 10% The same formulas, Eqs. ifT^ 
and H13(l , apply for the baryons with Sc replaced by St ■ If 
either the CDM or the baryon number was created after 
the curvaton decayed then there would be no isocurvature 
perturbation in that quantity . If both were created 
after the curvaton decayed there would be no isocurva- 
ture modes. 

There is no immediately compelling particle physics 
model for the curvaton scenario so we consider 

nine basic scenarios depending on whether the CDM and 
baryons are generated before, by, or after curvaton decay: 

If both the CDM and baryon number is created af- 
ter the curvaton decay then there is no isocurvature 
perturbation: 



1. 



(14) 



2. 



This scenario is consistent with the data and indis- 
tinguishable from an inflation model with negligible 
tensor component. 

If the CDM is created before the curvaton decays 
and the baryon number after the curvaton decays 
then from Eqs. ifT^ and ^ 



B = -3- 



Rh 



(15) 



This scenario is ruled out at high significance. 

3. If the baryon number is created before the curvaton 
decays and the CDM after the curvaton decays then 
from Eq. 

B = -3. (16) 
This scenario is ruled out at high significance. 

4. If the CDM is created by the curvaton decay and 
the baryon number after the curvaton decays then 
from Eqs. (O and © 



Solving for r gives 



1 



l + {Rb/R,)B/i- 



(17) 



(18) 



5. If the baryon number is created by the curvaton 
decay and the CDM after the curvaton decays then 
from Eq. ifT^ 



B^i 
Solving for r gives 



1 - r 



l + B/Z 



(19) 



(20) 



6. If the CDM and baryons are both created before 
the curvaton decays then from Eqs. © and (|12|l 



B 



3 

Rb 



(21) 



This scenario is ruled out at high significance. 

7. If CDM and baryons were both created by the cur- 
vaton then from Eqs. and (|13|l 



B 



Solving for r gives 



3 1 



Rh 



1 



l + RbB/?,' 



(22) 



(23) 



8. If the CDM is created before the curvaton decay 
and the baryons are created by the curvaton decay 
then from Eqs. ©, and (O 



B = 

Solving for r gives 



3(i?6 



rRb 



Rb 



l + RbB/5' 



(24) 



(25) 



9. If the CDM is created by the curvaton and the 
baryons are created before the curvaton decays we 
have from Eqs. I©, and (HSJ 



B = 



HRc - ' 
rRb 



Solving for r gives 



Rc 

l + RbB/3' 



(26) 



(27) 



For the cases that are not immediately ruled out we ob- 
tain a constraint on r. The posterior probability distri- 
bution for this quantity can easily be constructed from 
the Monte Carlo samples, and a plot of it is shown in 
Fig.|31for the various cases. The peaks at r = 1 are when 
there are no isocurvature modes. The curves which peak 
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FIG. 3: Plots of the un-normalized posterior probability dis- 
tribution for r ~ pcurvaton/ptotai when the curvaton decays. 
The distributions are for the numbered scenarios described 
in the text: (4) CDM created by curvaton decay and baryon 
number after curvaton decay (green long dashes), (5) baryon 
number created by curvaton decay and CDM after curvaton 
decay (cyan dash-dot line), (7) both CDM and baryon number 
created by curvaton decay (black solid line) , (8) CDM created 
before curvaton decay and baryon number by curvaton decay 
(blue dotted line) , (9) baryon number created before curvaton 
decay and CDM by curvaton decay (red short dashed line) . 



at r ~ i?h and r ^ Rc are when compensating baryon 
and CDM isocurvature modes are created before and by 
curvaton decay, giving a total effective isocurvature per- 
turbation close to zero. 

The amount of non-Gaussianity in the CMB is depen- 
dent on r with the conventional governing parameter 



assuming /„; ^ 1. Using this equation we can convert 
the likelihood plots for r into those for /„; as is shown in 
Fig. ^ The values near /„; ~ 1 should not be taken too 
seriously as there will be additional second order non- 
Gaussian contributions from fields other than the curva- 
ton. The current one year WMAP data has /„; < 134 
(95%) and is predicted to reach /„; < 80 (95%) with 
the four year WMAP data So if WMAP eventually 
detects non-Gaussianity it will rule out all the models 
considered here. The Planck satellite is predicted to ul- 
timately be able to detect /„; > 5 . If this is reaHzed 
Planck will be able to distinguish between the case where 
the CDM is created before curvaton decay and the baryon 
number by curvaton decay and the other possibilities. 



FIG. 4: Plots of the un-normalized posterior probability dis- 
tribution for the amount of non-Gaussianity, /„;. The line 
styles are the same as in Fig. |3 



IV. CONCLUSIONS 



The curvaton model provides a simple scenario that 
can give rise to correlated adiabatic and isocurvature 
modes of similar size. The current data do not favor 
a large isocurvature contribution, but a significant am- 
plitude is still allowed. 

We point out that the CDM and baryon isocurvature 
modes differ only by the addition of an observationally 
null mode in which the two perturbations compensate. 
The CDM isocurvature mode can therefore be treated 
as a scaled baryon isocurvature mode. A simple ana- 
lytical approximation for the effect of mixtures of large 
scale isocurvature and adiabatic perturbations on the 
CMB temperature anisotropy was given. Numerically, 
we found that the data was consistent at the two sigma 
level with the presence of an effective correlated baryon 
isocurvature perturbation of about 50% the magnitude of 
the adiabatic perturbation. The individual baryon and 
CDM isocurvature modes can be even larger if they com- 
pensate each other. Models in which either the baryon 
number or CDM was created before the curvaton dom- 
inated the energy density are ruled out unless counter- 
balanced by the other species being created by the cur- 
vaton decay. The levels of non-Gaussianity expected for 
the various scenarios were evaluated and in the case of 
the CDM being created before the curvaton decayed and 
the baryon number by the curvaton decay, could be high 
enough to detectable by the Planck satellite. 
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APPENDIX A: COVARIANT PERTURBATION 
EQUATIONS 

The covariant approach to cosmological perturbation 
theory gives a set of gauge invariant equations in which 
all the terms are covariant and have a physical interpre- 
tation [ssL l66j . The quantities can be calculated in any 
frame (labelled by a 4- velocity Ua) and the equations re- 
main the same. Individual quantities measuring a par- 
ticular perturbation do in general depend on what frame 
Ua is used to calculate them, so when talking about (for 
example) a density perturbation it is important to make 
clear what frame one is referring to. 

The spatial gradient of the 3-Ricci scalar '^^^TZ van- 
ishes in a homogeneous universe, and rja = ^SDa'^^^TZ is 
a natural covariant measure of the scalar curvature per- 
turbation in some frame with 4- velocity Ua- Here S is 
the scale factor and Da is the spatial covariant derivative 
orthogonal to Ua (we use the signature where Uau"" — 1). 
Other covariant quantities useful for studying perturba- 
tions are defined in 62, 66] , along with derivations of the 
equations of General Relativity that relate them. Here 
we only consider scalar modes at linear order in a spa- 
tially flat universe^, and perform a harmonic expansion 
as described in leaving the fc-dependence of scalar 
quantities implicit. For example we describe the curva- 
ture perturbation by the scalar harmonic coefficient rj. 

Frame invariant quantities can be constructed from 
combinations of covariant quantities that depend on the 
choice of frame Ua- These often have an interpretation 
in terms of the value of a particular quantity in some 
specified frame. In particular 



^ 1 



(Al) 



where a is the scalar shear and H = SSV^Ua is the con- 
formal Hubble parameter, is proportional to the curva- 
ture perturbation in the zero shear frame (the Newtonian 
gauge). The acceleration A in the zero shear frame 



* = -A + {a' + na)/k 



(A2) 



The equations given here generahze trivially to a non-flat uni- 
verse by the substitution — » »7/(l — 3K/k^). 



defines a second frame invariant quantity, which is related 
to $ by 



$ + ^ = - 



(A3) 



where H is the anisotropic stress. The Weyl tensor is the 
part of the Riemann tensor which is not determined by 
the local stress-energy, and defines a frame independent 
scalar potential^ cj) |6y| which is related to the above via 



= i(*-$). 

We define a frame invariant curvature perturbation 

7] HSp 



c 



p' 



(A4) 



(A5) 



proportional to the curvature perturbation in the uniform 
density frame. Here 5p = 5p^ = ^ ■ p^A^ is the total 
density perturbation. This is related to the comoving 
curvature perturbation 



X ■ 



Hq 



Hp + p) 



2 7^-1$' - 1- 

o — 

3 1 + w 



(A6) 



by X = — C — "HSp/ p' , where 5p is the comoving density 
perturbation and q ~ '^^{pi +Pi)vi is the total heat flux 
and Vi are the velocities. The Poisson equation relates 
the density and potential via fc^<I> = ^nS^Sp. It follows 
that for adiabatic modes where x is non-zero initially 
X ~ — C on large scales. 

A local scale factor S can be defined (up to an initial 
value) by integrating the local expansion rate 7i, and 
the quantity ha = DaS = SDaS/S (scalar harmonic 
coefficient h) provides a measure of the perturbation to 
local volume elements. The derivative h' with respect to 
conformal time r is unambiguously defined, and describes 
the rate of change of local volume element perturbations. 
In the frame in which h' is zero fractional perturbations in 
number densities of conserved species remain constant if 
there are no matter flows. The evolution of the curvature 
perturbation is given by 



1] 



2h' 



-k<j 



(A7) 



so on large scales the h' — frame coincides with the 
77' = frame. Thus 77 is conserved on large scales in 
the frame in which number density perturbation fractions 
are constant 0| . This result is purely a result of linear 
torsionless spacetime geometry. 

The time evolution of the local scale factor pertur- 
bation sources growth of density perturbations of unin- 
teracting conserved species via the energy conservation 
equation 

5p'^ + m{5p, + 5p,) + k{p,+p,)v, - -ih'{p,+pi), (A8) 



The <I> of Ref. 1661 has a different sign convention where <t> ; 



8 



where Spi is the pressure perturbation. The h' = frame 
therefore coincides with the 6p^ — frame on large scales 
if 5pi = in the Spj^ = frame. For a particular species 
one can define the curvature perturbation in the frame 
in which its density is unperturbed 



where in the absence of energy transfer p'^ = —i'H{p 
The evolution equation that follows from Eqs 

(Ixtii is m 



(A9) 



and 



n 



Pi +pi 



Spi - —5p, 



P 



m 

3 ' 



(AlO) 



where Vi = + cr is the Newtonian gauge velocity. If 
there is an equation of state pi = Pi{pi) the first term 
on the right hand side is zero, and the Ci are therefore 
constant on large scales where kr <^1. If the equation of 
state parameter Wi = Pi/ Pi is constant this implies that 
the fractional density perturbations in the unperturbed 
curvature frame evolve as 



A' 



(All) 



and hence the are also conserved on large scales. The 
curvature perturbation in the frame in which the total 
energy is unperturbed is given from the by Eq. |(lj. In 
the frame in which the acceleration A — (and hence Ua 
coincides with the CDM velocity) 77^ = — ??/2, h'^ = 6h' 
where rjs and h'^ are the synchronous gauge quantities 
(e.g. see [67|). 



APPENDIX B: ISOCURVATURE INITIAL 
CONDITIONS 

In the early radiation dominated era there are in 
general five regular solutions to the perturbation equa- 
tions assuming there is only one distinct species of 
cold dark matter. If there are several species of dark 
matter the additional modes are unobservable without 
measuring the distinct dark matter species directly. Per- 



forming a series expansion in conformal time r, the Fried- 
man equation gives 



n, 



,^22 

UJT H UJ T 

4 



0[Kujt'^) 



(Bl) 



where uj = D,„i'H{)/ ^JT^'+Tl^ with TLq the Hubble pa- 
rameter today and 17,; the density today in units of the 
critical density. At lowest order in the tight coupling 
expansion, assuming the baryons and dark matter have 
negligible pressure, the CDM isocurvature mode at early 
times is 



A, 
A, 



1 



1 RciAR^ - 15)cjfc2r3 



72 



2i?^ + 15 



- _ 4 . _ 5 RcUjPt^ 

3 6 2(i?, + 15) 

1 Rc{4:R^ - Idjujkr'^ 
24 2R^ + 15 

15 Rrujkr'^ 



= V, = Vb = — 



2R^ + 15 
1 i?c(4i?^ - 15)t^T 
8 2R^ + 15 
1 i?c(4i?^ + l5)ujT 
8 2i?^ + 15 



(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 



where equalities apply at the given order in t. The 
baryon isocurvature mode is given by subtracting the 
observationally null mode Ac = —Rb^^b/Rc from the 
above solution. Series solutions for the adiabatic and 
isocurvature modes to any order are easily computed 
using computer algebra packages, for a Maple deriva- 
tion of the solutions in the zero acceleration frame see 
|http : / / Ceunb . info/theory .html The above solution 
was calculated by constructing the frame invariant quan- 
tities above from the quantities in the zero acceleration 
frame. 

The Ai are constant to order (kr)'^. However the low- 
est order terms in the velocities are of order (a;T)(/cT), 
demonstrating explicitly that the assumption that Vi = 
const X kr + 0{{kT)^) in is incorrect for isocurvature 
modes. 



[1] G. Hinshaw et al. (2003), astro-ph/0302217. 
[2] H. V. Peiris et al. (2003), astro-ph/0302225. 
[3] A. Liddle and D. Lyth, Cosmological Inflation And Large- 

Scale Structure (Cambridge University Press, 2000). 
[4] L. A. Kofman and A. D. Linde, NucL Phys. B 282, 555 

(1987). 

[5] D. Polarski and A. A. Starobinsky, Phys. Rev. D 50, 6123 

(1994), astro-ph/9404061. 
[6] J. Garcia-Bellido and D. Wands, Phys. Rev. D 53, 5437 

(1996), astro-ph/9511029. 
[7] A. Linde and V. Mukhanov, Phys. Rev. D 56, 535 (1997), 

astro-ph/9610219. 



[8] M. Sasaki and T. Tanaka, Prog. Theor. Phys. 99, 763 

(1998), gr-qc/9801017. 
[9] D. Langlois, Phys. Rev. D 59, 123512 (1999), astro- 

ph/9906080. 

[10] C. Gordon, D. Wands, B. A. Bassett, and R. Maartens, 

Phys. Rev. D 63, 023506 (2001), astro-ph/0009131. 
[11] J.-C. Hwang and H. Noh, Phys. Lett. B 495, 277 (2000), 

astro-ph/0009268. 
[12] D. Wands, N. Bartolo, S. Matarrese, and A. Riotto, Phys. 

Rev. D 66, 043520 (2002), astro-ph/0205253. 
[13] S. Tsujikawa, D. Parkinson, and B. A. Bassett (2002), 

astro-ph/0210322. 



9 



A. A. Starobinsky, S. Tsujikawa, and J. Yokoyama, Nucl. 
Phys. B610, 383 (2001), astro-ph/0107555. [41 
N. Bartolo, S. Matarrese, and A. Riotto, Phys. Rev. D64, 
123504 (2001), astro-ph/0107502. [42 
S. Groot Nibbelink and B. ,]. W. van Tent, Class. Quant. 
Grav. 19, 613 (2002), hep-pli/0107272. [43 
S. Mollerach, Phys. Rev. D 42, 313 (1990). 
D. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002), [44 
hep-ph/0110002. 

D. H. Lyth, C. Ungarelh, and D.Wands, Phys. Rev. D [45 
67, 023503 (2003), astro-ph/0208055. 
T. Moroi and T. Takahashi, Phys. Lett. B 522, 215 [46 

(2001) , hep-ph/0110096. [47 
T. Moroi and T. Takahashi, Phys. Rev. D 66, 063501 [48 

(2002) , hep-ph/0206026. [49 
T. Moroi and H. Murayama, Phys. Lett. B 553, 126 

(2003) , hep-ph/0211019. [50 
K. Dimopoulos and D. Lyth (2002), hep-ph/0209180. 

M. Bastero-Gil, V. Di Clemente, and S. F. King (2002), [51 
hep-ph/0211011. 

M. Postma, Phys. Rev. D 67, 063518 (2003), hep- [52 
ph/0212005. 

K. Enqvist, S. Kasuya, and A. Mazumdar, Phys. Rev. [53 
Lett 90, 091302 (2003), hep-ph/0211147. [54 
K. Dimopoulos, G. Lazarides, D. Lyth, and R. Ruiz de 
Austri (2003), hep-ph/0303154. [55 
K. Dimopoulos, D. H. Lyth, A. Notari, and A. Riotto [56 
(2003), hep-ph/0304050. 

K. Enqvist and M. S. Sloth, Nucl. Phys. B626, 395 [57 
(2002), hep-ph/0109214. 

M. S. Sloth, Nucl. Phys. B 656, 239 (2003), hep- [58 
ph/0208241. 

V. Bozza, M. Gasperini, M. Giovannini, and [59 
G. Veneziano, Phys. Lett. B543, 14 (2002), hep- 
ph/0206131. [60 

E. J. Copeland, R. Easther, and D. Wands, Phys. Rev. 
D56, 874 (1997), hep-th/9701082. [61 
J. E. Lidsey, D. Wands, and E. J. Copeland, Phys. Rept. 

337, 343 (2000), hep-th/9909061. [62 
G. Efstathiou and J. R. Bond, Mon. Not. R. Astron. Soc. 
218, 103 (1986). [63 
R. Stompor et al., Astrophys. J. 463 (1996). [64 
E. Pierpaoli, J. Garcia-Bellido, and S. Borgani, JHEP 
10, 015 (1999), hep-ph/9909420. [65 
M. Kawasaki and F. Takahashi, Phys. Lett. B516, 388 
(2001), hep-ph/0105134. [66 
K. Enqvist, H. Kurki-Suonio, and J. Valiviita, Phys. Rev. 
D62, 103003 (2000), astro-ph/0006429. [67 
D. Langlois and A. Riazuelo, Phys. Rev. D 62, 043504 
(2000), astro-ph/9912497. 

M. Bucher, K. Moodley, and N. Turok, Phys. Rev. D 62, 



083508 (2001), astro-ph/9904231. 

R. Trotta, A. Riazuelo, and R. Durrer, Phys. Rev. Lett. 
87, 231301 (2001), astro-ph/0104017. 
R. Trotta, A. Riazuelo, and R. Durrer, Phys. Rev. D 67, 
063520 (2003), astro-ph/0211600. 

M. Bucher, K. Moodley, and N. Turok, Phys. Rev. Lett. 

87, 191301 (2001), astro-ph/0012141. 

L. Amendola, C. Gordon, D. Wands, and M. Sasaki, 

Phys. Rev. Lett. 88, 211302 (2002), astro-ph/0107089. 

N. Bartolo and A. R. Liddle, Phys. Rev. D65, 121301 

(2002), astro-ph/0203076. 

L. Verde et al. (2003), astro-ph/0302218. 

A. Kogut et al. (2003), astro-ph/0302213. 

C. L. Kuo et al. (2002), astro-ph/0212289. 

W. Percival et al., MNRAS 337, 1068 (2002), astro- 
ph/0206256. 

W. L. Freedman et al., Astrophys. J. 553, 47 (2001), 
astro-ph/0012376. 

S. Buries, K. M. NoUett, and M. S. Turner, Astrophys. 

J. 552, LI (2001), astro-ph/0010171. 

A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002), 

astro-ph/0205436. 

J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

J. M. Bardeen, P. J. Steinhardt, and M. S. Turner, Phys. 

Rev. D 28, 679 (1983). 

D. H. Lyth, Phys. Rev. D 31, 1792 (1985). 

H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 
78, 1 (1984). 

V. Mukhanov, H. Feldman, and R. Brandenberger, Phys. 
Rep. 215, 203 (1992). 

D. Wands, K. A. Malik, D. H. Lyth, and A. R. Liddle, 
Phys. Rev. D 62, 043527 (2000), astro-ph/0003278. 

K. H. Malik, D. Wands, and C. Ungarelh, Phys. Rev. D 
67, 063516 (2003), astro-ph/0211602. 
U. Seljak and M. Zaldarriaga, Astrophys. J. 469, 437 
(1996), astro-ph/9603033. 

A. Lewis, A. Challinor, and A. Lasenby, Astrophys. J. 
538, 473 (2000), http://camb.info astro-ph/9911177. 
A. Lewis and A. Challinor, Phys. Rev. D 66, 023531 
(2002), astro-ph/0203507. 

E. Komatsu et. al (2003), astro-ph/0302223. 

E. Komatsu and D. N. Spergel, Phys. Rev. D 63, 063002 
(2001), astro-ph/0005036. 

G. F. R. Ellis, D. R. Matravers, and R. Treciokas, Ann. 
Phys. 150, 455 (1983). 

A. Challinor and A. Lasenby, Astrophys. J. 513, 1 (1999), 
astro-ph/9804301. 

C.-P. Ma and E. Bertschinger, Astrophys. J. 455, 7 
(1995), astro-ph/9506072. 



